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ABSTRACT 

In this paper we derive ages and masses for 276 clusters in the merger galaxy 
NGC 3256. This was achieved by taking accurate photometry in four wavebands from 
archival HST images. Photometric measurements are compared to synthetic stellar 
population (SSP) models to find the most probable age, mass and extinction. The 
cluster population of NGC 3256 reveals an increase in the star formation rate over the 
last 100 million years and the initial cluster mass function (ICMF) is best described 
by a power law relation with slope a = 1.85 ± 0.12. 

Using the observed cluster population for NGC 3256 we calculate the implied 
mass of clusters younger than 10 million years old, and convert this to a cluster 
formation rate over the last 10 million years. Comparison of this value with the star 
formation rate (SFR) indicates the fraction of stars found within bound clusters after 
the embedded phase of cluster formation, T, is 22.9%±gJ for NGC 3256. We carried 
out an in-depth analysis into the errors associated with such calculations showing that 
errors introduced by the SSP fitting must be taken into account and an unconstrained 
metallicity adds to these uncertainties. Observational biases should also be considered. 

Using published cluster population data sets we calculate T for six other galaxies 
and examine how T varies with environment. We show that T increases with the star 
formation rate density and can be described as a power law type relation of the form 
T(%) = (29.0 ± 6.O)E° j, 4 ± o o4 (M yr- 1 kpc~ 2 ). 

Key words: galaxies: structure - galaxies: stellar content - stars: formation - stars: 
clusters 



1 INTRODUCTION 

Galaxy mergers produce some of the most extreme star for- 
mation events in the universe, inducing starbursts with in- 
tensities (i.e. star-formation rates; SFRs) orders of magni- 
tude larger than in quiescent galaxies. In the distant universe 
these can be seen as Hyper-luminous Infrared galaxies (e.g. 
|Verma et al.|2002| with SFRs exceeding a few thousand solar 
masses per year. In the local universe, ongoing mergers like 
Arp 220 e.g. (e. g|Scoville et al.|1998||Wilson et al.|2006| > and 
NGC 3256 (e.g. |Zepf et al.|1999||Trancho et al.|2007b| > have 
SFRs of a few tens to hundreds of solar masses per year, 
and due to their proximity, can be resolved with Hubble 
Space Telescope imaging. In these systems large numbers of 
super-star clusters are often found , with ages of a 1-500 Myr 
and masses up to ~ 8 x 1O 7 M0 ( Maraston et al.|2004 |. We 
can think of star clusters as simple stellar populations which 
can be modelled relatively easily. This, combined with their 
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high surface brightness that allows them to be easily de- 
tected, and their long lifetimes, make them ideal tracers of 
the star formation history of a galaxy. 

NGC 3256 is a relatively nearby starburst galaxy that 
is clearly the result of a recent galactic merger. It displays 
two prominent tidal tails, thought to have been produced 
during the first encounter between two spiral galaxies ap- 
proximately 500 million years ago ( Zepf et al.|1999 English 
et al.||2003[ ). |Zepf et a!| ( |1999[ ) catalogued over 1000 young 
star clusters in the main body of the merger using HST 
imaging. Trancho et al. (2007b I discovered three massive 
(1 — 3 x 10 a Mq) clusters within one of the tidal tails, whose 
ages and velocities places their formation within the tidal de- 



bris. Additionally, Trancho et al. (2007a I (hereafter T07a) 



studied a further sample of 23 clusters in the main body of 
the remnant and found that the clusters had metallicities 
of ~ 1.5 Zq, masses in the range 2 — 4 x 1O 5 M0 and ages 
from a few to 150 Myr. The current SFR of NGC 3256 is 
~ 50 MQyr -1 ( Bastian 20081, however it's SFR appears to 
have been increasing for the past ~ 200 Myr (see § 2.5 1 and 



2 Q. E. Goddard, N. Bastian, and R. C. Kennicutt 



Trancho et al. (2007a I argue that it is likely to continue to 



increase in the future. 

The age distribution of clusters reveals the underlying 
star formation history of the host galaxy. The mass distribu- 
tion of clusters gives information on the cluster initial mass 
function (CIMF) which is often approximated by a power- 
law of the form Ndm oc M~ a dm, with a ~ 2 (Elmegreen 



fc Efremov|[T997l |Zhang fc FaHp999l |de Grijs et aTf 2003c 



Bik et al.||2003 1. More recently, it has been shown that the 
CIMF may have a truncation at high masses, being well de- 
scribed by a Schechter function ( |Schechter|1976[ ) which is a 
power-law in the low-mass regime and has an exponential 



cutoff above a given mass, M*(Gieles et al. 2006a Bastian 
[2008} |Larsen|2009| |GleTes|2009b| ). 

In order to obtain accurate cluster ages and masses ei- 
ther spectroscopy or photometry may be used. Although 
spectra yield more accurate results, photometry is more ap- 
plicable to large sample sets. In order to break the age- 
dust degeneracy, a cluster must be observed in at least four 



broad photometric bands and cover the Balmer break (An- 



ders et al. 2004a I. By comparing the observed colours and 



magnitudes of each cluster to synthetic stellar population 
(SSP) models, we can estimate the age, mass and extinction 
of each cluster. This technique has been used on several clus- 
ter populations, the Antennae ( Fall et al.|2005 Anders et al 



2007), NGC 1569 (Anders et al 



2005b) 



2004a I & M51 (Bastian et al 



If all stars are formed in clusters, then only a small per- 
centage of them are still within clusters at the end of the 



embedded phase (Lada & Lada 20031. For up to the first 



three million years of a cluster's lifetime it remains embed- 
ded in the progenitor molecular cloud until stellar winds, 
ionising flux from massive stars, and stellar feedback expel 



the remaining gas (see Goodwin]|2009 for a recent review). 
In expelling this gas many clusters become unbound and 



subsequently disperse into the surrounding medium (Lada 
fc Lada|2003"l ), a process commonly described as 'infant mor- 



tality'. Even if clusters do survive the embedded phase in- 
tact other disruption mechanisms may come into effect, such 
as stellar evolutionary mass loss, two body relaxation and 
GMC encounters. These mechanisms generally operate over 
longer time-scales, being a few 10s of Myr for stellar evo- 
lution (e.g. iBastian et aLl|2009[ ) and 100s of Myr for GMC 
encounters ( Gieles et al. 2006c I, except when the GMC num- 
ber density is particularly high (e.g. Lamers et al.|2005 l. 



After a cluster disrupts the stars disperse and become 
part of the galactic background. Comparing the fraction 
light emitted from clusters to the total (i.e. background plus 
clusters) galactic emission gives an indication of the fraction 
of stars within bound clusters. Meurer et al. ( 1995 1 esti- 
mated that 20 — 50% of UV light in starburst galaxies comes 



from star clusters. Zepf et al. ( 1999 1 calculated the fraction 
of light from clusters in the B band to be 15 — 20% and 
half that in the I band for the galaxy NGC 3256. However 
the most comprehensive set of results comes from |Larsen| 
& Richtler ( 2000 1 , listing the fraction of light from clusters 



in both the U-band (T L (U)) and V-ban ds (T £ (V)) for 32 
galaxies with varying star formation rates. [Larson & Ri chtler| 
(20001 found that T_l(U) increases with the star formation 
rate of the host galaxy and a stronger correlation with the 
star formation rate surface density, indicating that the host 



environment may influence the mode of star and/or cluster 
formation and how likely a cluster is to survive. 

Measuring the fraction of light from clusters is a rela- 
tive easy calculation assuming foreground stars can be elim- 
inated from the sample. However, it is presently unknown, 
but calculable, how these values may be influenced by the 
presence of bursts of star formation in the past and differ- 
ential extinction of younger clusters. Hence the fraction of 
light observed in clusters will be a (possibly complicated) 
combination of the fraction of stars formed in clusters, the 
star-formation history of the galaxy, the cluster disruption 
time-scale/rate, and the difference in the amount of extinc- 
tion towards clusters and the field. 

In this paper we attempt to improve the situation by de- 
riving ages and masses for clusters in the galaxy NGC 3256, 
which in turn is used to calculate a cluster formation rate 
(CFR) over the last ten million years. Comparing this in- 
ferred CFR to the star formation rate (SFR, measured by 
Ha fluxes or infrared luminosities) we compute the fraction 
of stars in clusters younger than ten million years which 
have survived the embedded phase intact, a value hereafter 



referred to as T (Bastian 20081. We pay particular atten- 



tion to the possible sources of uncertainty which may affect 
these calculations. Using other data sets of cluster ages and 
masses for different galaxies we perform the same calcula- 
tion for an additional 6 galaxies in an attempt to find how 
r varies with environment. 

Throughout this paper we will define a "cluster" as a 
gravitationally bound, centrally concentrated group of stars 
that can be identified on high resolution optical imaging. 
Hence, clusters refer to objects that have survived the transi- 
tion from being embedded in their natal GMC to an exposed 
state. We assume that this process happens approximately 
3 Myr after the cluster forms. We adopt a Hubble constant 
of H = 70 km s" 1 Mpc -1 , which places NGC 3256 at a dis- 
tance of 36.1 Mpc given it has a recession velocity relative 
to the local group of +2804 ± 6 km s~\ This corresponds 
to a distance modulus of 32.79. 

This paper is structured as follows, we begin in § [2] 
by introducing the dataset for NGC 3256 and the methods 
we used to determine the age and masses of clusters in the 
galaxy. In §[3] we detail how we calculated F for NGC 3256 
and go on to carefully examine all the possible sources of 
error in making such a calculation. In §[4] we introduce other 
datasets taken from the literature and calculate F for each. 
In § [5] we search for trends with T and galactic properties 
and discuss the implication. Finally, in § [6] we present our 
conclusions and summarise our main results. 



2 OUR STUDY OF NGC 3256 
2.1 The Data 

We used archival Hubble Space Telescope (HST) Advanced 
Camera for Surveys (ACS) images of NGC 3256 across four 
filters, F330W (PI: Holland Ford; ID: 9300), F435W (PI: 
Alex Fillippenko; ID: 10272), F555W (ID: 9300) and F814W 
(ID: 9300). Details of all the images used can be found in Ta- 
ble [l] The High Resolution Camera (HRC) images all cover 
the same area apart from the HRC-F^^H^ image which is 
offset. To ensure we were examining the same area in all four 
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Table 1. Details of images 



Filter 


Camera 


Exposure 


Aperture 


Name" 






Time (s) 


Correction 




F330W 


HRC 


2000 


0.418 


U 


F435W 


HRC 


840 






F555W 


HRC 


760 


0.397 


V 


F814W 


HRC 


660 


0.405 


1 


F435W 


WFC 


1320 


0.328 


B 



a We refer to the filters by their closest Cousins- Johnson filter 
name, however no transformations have been applied. 

bands we used the ACS Wide Field Camera (WFC) F435W 
(PI: Aaron S. Evans; ID: 10592) image over the HRC F435W 
image. The WFC image has the advantage of a longer ex- 
posure time, though the pixel scale is slightly larger, 0.05" 
compared to 0.027". Throughout this paper we will refer to 
the F330W, F435W, F555W, and F814W in the Cousins- 
Johnson U, B, V and I notation, although we stress that no 
transformations have been applied. 



2.2 Cluster Selection & Photometry 

We selected star clusters using the SExtractor program 
( |Bertin fc Arnouts|p96| , run on the B-band WFC image 
in order to maximise the number of clusters detected. We 
changed the SExtractor detection parameters to maximise 
the number of objects found, attempting to avoid blending 
of objects in the most crowded regions. We set a minimum 
detection limit as 4cr above the local background as deter- 
mined by SExtractor , and a minimum object area of 8 pix- 
els. Co-ordinates of the clusters were translated from the B 
band image into co-ordinates for the other bands using the 
GEOMAP and GEOXYTRAN routines within IRAF. 

Photometry across all images was carried out using the 
A PER routine from the DAOPHOT package in IRAF. We 
did not use the photometry from SExtractor to ensure a 
consistent and unbiased method in all wavebands. We used 
an aperture of radius 7.0 pixels for the HRC images and 3.5 
pixels for the WFC image. An annulus with radius 8 pix- 
els and a thickness of one pixel was used to measure the 
local background in the HRC images, radii and thickness 
of the annulus were adjusted accordingly for the WFC im- 
age. Aperture corrections were calculated based on several 
bright, isolated and spatially resolved clusters, by compar- 
ing their flux in an aperture of 30 (HRC) pixels relative to 
our adopted 7 (HRC) pixel aperture. The resulting aperture 
corrections are shown in Tabled] 

We corrected the photometry for foreground galactic 
extinction, using the extinction law of Savage & Mathis 
(19791, and a dust extinction value of Ay = 0.122, taken 



from Schlegel et al. (19981 



In total we measured fluxes for 904 clusters in NGC 
3256 across the U, B, V and I bands. Figure [I] shows a 
colour-colour diagram of clusters within a magnitude limit 
(m_g < 21.0). Figure [5] shows a colour-composite image of 
NGC 3256. To avoid overcrowding we have only marked the 
positions of the 276 cluster that have good fits to SSP mod- 
els, more details on this selection can be found in section 
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Figure 1. Colour-Colour diagram for clusters below the magni- 
tude limit mg < 21.5 (simply to avoid overcrowding the figure). 
Clusters have not been reddening corrected. We have only in- 
cluded clusters that met our criteria for a good fit to the SSP 
models, as explained in section [2. 3| The lines show solar metallic- 
ity SSP models, the solid line is based on the Padova isochrones 
whilst the dashed line is based on the Geneva isochrones. The 
red dashed segment shows the 'red-loop', few clusters are ever 
associated with this feature. 
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Figure 2. Colour composite image of of NGC 3256. We have 
marked on the positions of the 276 clusters classified as good fits 
in section [2. 31 



2.3 Determination of Cluster Parameters 

In order to determine the age, mass and local extinction 
for each cluster the photometric values were compared to 
those of cluster evolution models, in a similar fashion to 
|Bik et al.| ( |2003[ ) and |Bastian et ah] (2005a} . We used the 
most recent GALEV simple stellar population (SSP) models 
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( Anders fc Fritze-v. Alvensleben|2003| | Anders et al.|2004a |. 
Although SSP models exist based on both the Padova and 
Geneva isochrones, we have only used the Padova tracks to 
determine cluster parameters. The Geneva tracks display a 
feature known as the red loop and tend to produce a poor 



fit to data compared to the Padova tracks (Whitmore & 



Zhang 2002). The GALEV models have the advantage of 



being produced with colours matching the HST niters, this 
avoids errors associated with converting between HST fil- 
ters and the standard Cousins-Johnson system. In addition, 
the GALEV models include gaseous emission lines and con- 
tinuum emission. The youngest data points in the GALEV 
models are 4Myr and 8Myr, this can lead to poor fitting of 



the youngest clusters ( jBastian et al. 2005a I. To avoid this we 
linearly interpolate between the youngest model ages for all 
filters. This produces SSP tracks which are regularly sam- 
pled at intervals of 0.02 in log time. The determined age of 
a young cluster is not bound to 4 or 8Myr in this case and 
can take a more precise value with a better fit over the four 
wavebands. 

We have used a three-dimensional maximum likehood 



fitting method developed by Bik et al. (20031. The method 



is described in detail there, and so we only summarise the 
process in this paper. We also note that this method has 
been tested against colour-colour methods and has shown 
itself to be superior ( de Grijs et al.|2003b Parmentier et al. 
2003 ). 

In simple terms the method works as follows: the 
GALEV SSP models give a grid of colours as a function 
of age. We then apply extinction to each model in steps of 
0.02 in E(B-V) to extend the grid. Each cluster then has its 
observed spectral energy distribution compared to the grid 
using a minimum % 2 test. The model with the age and ex- 
tinction which produces the lowest % 2 ls chosen, and a range 
of acceptable values is calculated. As the GALEV model 
magnitudes are scaled to a mass of IO 6 M the difference 
between the observed magnitudes and model magnitudes 
can be converted to a mass, given the distance modulus and 
extinction. The accuracy of this fitting method is not stud- 



ied directly here as it has been reported previously (Bastian 
et al.||2005a I, for the majority of the sample errors on the 
age and mass are less than 0.3 dex. 

After obtaining the best fit for each cluster we sought 
to determine if the fit was good or not. In fitting cluster 
parameters a x 2 value is calculated and this is one way of 
identifying a good fit. However for bright objects with low 
photometric errors the \ 2 value may be high and thus result 
in a bad fit. To avoid such instances we use the standard 
deviation, defined in Equation [T] 



nfilte 



(1) 



The a value is summed over all four wavebands. The 
spectral energy distributions were compared to the values 
predicted by the best fit models and compared by eye. This 
was done to establish a minimum value for a 2 , a^ in . We 
also computed a histograms of a 2 for all sources and ap- 
proximating the distribution with a gaussian and took the 
\ point. This value was very close to an appropriate value 
for <7m in determined by visual inspection at 0.03. 

In the resulting analysis of the properties of cluster sys- 



tem of NGC 3256 we only consider those clusters which pass 
the following criteria: 

(i) detected in all four wavebands, at least 5a above the 
background; 

(ii) uncertainties in the magnitude of all wavebands ^0.2 
mag; 

(iii) is well fitted by an SSP model {a sC 0.03). 

The number of clusters which passed these rigourous se- 
lection criteria depends on the metallicity of the SSP model 
used. With a solar metallicity Padova model we accept 276 
objects, with a twice solar metallicity we accept 285 objects. 
Details of the 276 clusters accepted using the solar metal- 
licity Padova tracks can found with the electronic version 
of this article. It should be noted that the higher number 
of accepted objects produced by the twice solar metallic- 
ity models is not an indication of the global metallicity of 
the galaxy. Metallicity is very hard to fit accurately with 



only photometry, as shown by Bastian et al. (2005a I. Tran- 
cho et al. (2007a I found the metallicity to range between 



LI - 1.7 Z for 23 clusters in NGC 3256. 



Trancho et al. (2007a I (hereafter T07a) presented op- 



tical spectroscopy of 23 clusters in NGC3256. We have at- 
tempted to compare our results with those of T07a, however 
there are several problems that we encountered. Firstly, the 
astrometry presented by T07a is not accurate enough to 
identify any one cluster from our study, we were forced to 
visually inspect images from the two studies and then at- 
tempted to find the most likely matching clusters. Secondly, 
several of the clusters presented by T07a are in fact com- 
plexes, comprised of several clusters. Out of the ten clusters 
in the T07a study that lie within the same field of view as 
our study we found that eight had corresponding clusters in 
our study which matched our criteria for a good cluster fit. 
The remaining two clusters from the T07a were both com- 
plexes located in extremely crowded regions of the galaxy 
and could correspond to any one of several clusters in our 
study. Of the eight matched clusters T07a states that six of 
these were emission line objects and have reported ages of 
less than 10 6 ' 8 years. We find that ages for seven of these 
clusters is also under 10 6 ' 8 years, the remaining object has 
an age between 10 6 ' 7 - 10 6 ' 9 which appears to be consistent 
with the findings of T07a. We have also be able to match 
two absorption line clusters from the T07a study, both of 
these clusters have reported ages that are consistent within 
their respective errors. 

2.4 Cluster Properties of NGC 3256 

We can determine the properties and history of NGC 3256 
by examining the clusters that reside within the galaxy. 
With the age, mass and extinction of the clusters known we 
can derive several important properties. Figure [3] shows the 
age/mass distribution for NGC 3256 derived from the solar 
metallicity Padova SSP models. Immediately obvious from 
Figure [3] is the lack of low mass, old clusters. This is due to 
the fact that clusters dim as they age and so eventually be- 
come fainter than our detection limit. We can also note the 
large number of clusters seen with ages below lOMyr and 
over a large range of masses. Although the cluster fitting 
method can create some observed structure in the age/mass 
diagram ( Gieles et al. |[2005 I it is unlikely to do so over all 
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Figure 3. The age-mass distribution for all clusters with a good 
fit, using solar metallicity Padova SSP models. The dashed lines 
represent mass cuts at 4.7 and 5.7 log(MQ), which are used in 
figure [5] and throughout the rest of the study. The solid black 
line represents the magnitude limit my = 23.5 based on the solar 
metallicity Geneva SSP models at the distance of NGC 3256 



masses at low ages. We can conclude that there is a genuine 
over-density of clusters with ages below lOMyr. This can be 
interpreted as a lack of older clusters which could be due in 
part to the disruption of clusters, a point we shall address 
when we look as the cluster formation history in § |2.5| 

Various mass limits are shown on Figure [3] as horizontal 
dashed lines. Due to the effects of old clusters becoming 
fainter, higher mass limits result in a population which is 
complete up to a larger range of cluster ages. For example a 
mass limit of log(M /Mq) > 4.5 would only be complete for 
clusters younger than lOMyr. A mass limit of log(M / Mq) > 
5.5 would however be complete for much older clusters, those 
less than ~ 200 Myr. 

We can compute the slope of the cluster mass function 
for clusters younger than 10 million years which we have 
shown in Figure [4] This assumes the number of clusters of 
a given mass follows a power law distribution, like that in 
Equation [2] n(M) is the number of clusters with mass M, 
X is a constant of normalisation and a is the slope of the 
power law. 



N{m)dm = \M a dm 



(2) 



When calculating the slope of a mass distribution for 
any type of object th e way in which the data are bi nned 
can affect the result (Maiz Apellaniz & Ubeda 20051. To 



avoid such biases we have used bins with an equal number of 
objects in each. The minimum and maximum limits to each 
bin then depend on the masses of objects within. We have 
chosen to calculate the mass function for young clusters with 
ages less than lOMyr, in order to avoid any mass dependent 
cluster disruption effects, resulting in a complete sample. 

This calculation produces the mass function shown in 
Figure [4] where the data points are shown in black and the 
line of best fit is in red. The slope is fit to data points with 
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Figure 4. Mass function for clusters younger than 10 million 
years old in NGC 3256, produced using bins with a constant num- 
ber of clusters per mass bin. The red line is a line of best fit for 
points beyond the low mass turn over. 



masses above the limit M/Mq > Zog(4.3), below which we 
see a turn over in the mass function. This common turn 
over in cluster mass functions is due to incompleteness in 
the sample for low mass clusters. The slope of the mass 
function for NGC 3256 was measured to be a = 1.85 ±0.12, 
which is similar to other observed cluster populations that 
are generally in the range of a = 1.8 — 2.2 (de Grijs et al. 
2003a| |McCrady fc Graham||2007 1. 



2.5 Cluster Formation History of NGC 3256 

Ages determined by fitting photometric observations to SSP 
models can lead to large uncertainties, especially for faint ob- 
jects with large photometric errors. In constructing a cluster 
formation history for NGC 3256 it is important to take the 
uncertainty in the cluster ages into account. 

We have constructed a 4?r distribution using the same 



method described by Gieles et al. (20071. We do not reit- 



erate the details laid out in Gieles et al. (20071, only the 
central ideas. Each cluster has a contribution to the over- 
all age distribution as defined by an asymmetrical Gaussian 
spread in log t. This is calculated based on the minimum 
and maximum allowed ages of each cluster, calculated by 
the SSP fitting routine. 

In Figure [5] we show the smoothed age distribution re- 



sults for two different lower mass limits 

-,5.7 



10 4 7 M Q (left 



panel) and 10 Mq (right). The red shaded area shows 
the la Poisson errors, estimated by counting the number of 
clusters in bins of width 0.25, corresponding to the mean 
uncertainty in the log age values of the clusters. 

Using two different lower mass limits in Figure [5] our 
sample of clusters is complete up to two ages, the 1O 4T M 
limit shown in the left panel is complete for clusters younger 
than 20 million years. The higher mass limit of 10 5 - 7 M 
shown on the right is complete for clusters younger than 
200 million years, both of these ages are shown as a vertical 
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Age (x) / yr Age (x) / yr 

Figure 5. Apparent rate of cluster formation as a function of time for NGC 3256 and 2 different lower cluster mass limits. Left panel 
uses a limit of 10 4 ' 7 Mq which is complete for clusters younger than 20 million years, represented by the solid vertical line. The right 
panel shows the same plot for clusters more massive than 10 5 ' 7 Mq, which is complete up to 200 million years. The red shaded area 
denotes the Icr Poisson errors as described in the text. 



line on each plot. In both panels the cluster formation rate 
declines for clusters older than the completeness limit, due 
in part to older clusters fading and becoming undetectable. 

Focussing on the high mass cut sample, we see that 
the cluster formation rate (CFR) has increased by a factor 
of ~ 10 during the past ~ 200 Myr. This increase is ex- 
pected due to the ongoing galactic merger which is presum- 



ably the cause of the ongoing starburst (Bastian et al. 2009 1 



This interpretation is also supported by the present SFR of 
NGC 3256 which is ~ 46 Moyr -1 , which is approximately 
ten times higher than the sum of the two progenitor spiral 
galaxies. As in the Antennae ( |Bastian et al.|20l )9;), we do not 
see evidence for a high degree of long duration (> 10 Myr) 
mass independent cluster disruption, as proposed by |Fall| 
|et al.| ( p005] ). However, our high mass cuts means that we 
are insensitive to any disruption in the lower mass end of 
the cluster mass function, which may be strongly effected 
by cluster disruption during galaxy mergers (Kruijssen et 
al. 2010 in prep.). 



3 CALCULATING THE RATIO OF STARS 
FORMING WITHIN CLUSTERS, T 

3.1 Initial Estimate 

Assuming we were able to detect all the clusters within a 
galaxy we would be able to establish the total mass of clus- 
ters which had recently formed. Comparing this mass with 
a measure of the star formation rate (SFR) over the same 



time scale gives the ratio of stars forming within clusters, 
hereafter referred to as V (Bastian 2008). 

In principle this calculation is simple but is made more 
complex by the incompleteness of our sample of clusters. 
We are not able to detect low mass clusters, even those with 
young ages. Any low mass cluster is also likely to have large 
photometric errors and thus be rejected as a poor fit. We 
know, however, our sample is likely to be complete for bright, 
massive, young clusters. We can then extrapolate the mass 
found in these clusters to calculate the total mass of all 
clusters. 

In order to calculate the ratio of mass found above a cer- 
tain mass limit we made sample cluster populations based on 
a power law distribution (stochastically sampled), like that 
in Equation [2] We made model populations with 500,000 
clusters with a lower mass limit of IOOMq, and an upper 
limit of 1O 11 M0. This upper mass limit may appear very 
high but with a power law distribution with index of —2.0 
(i.e. q = 2.0) we only expect 50 clusters above 10 6 Mq based 
on this upper mass limit. To avoid erroneous results we ig- 
nore clusters with masses larger than 1O 7 M0. The upper 
mass limit allows the model populations to be fully popu- 
lated for high masses. We also calculate an error associated 
with the fraction of mass found above a certain mass limit, 
this is the standard deviation of results from 200 model 
cluster populations. The error does not decrease with an 
increased number of model populations as it is limited by 
stochastic variations in the masses of high mass clusters, 
which are rare in each model population. We have used 
a power law index of —2.0 in our claculations although it 
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should be noted that an index of —1.85 was measured for 
NGC 3256. We assume an index of —2.0 as this is the value 
found for large samples of clusters, we do accept that this 
value may vary between galaxies and the implications of this 
are discussed further in § |3.2| 

At this point we note that sampling an ICMF to gener- 
ate a cluster population simply to calculate the mass above 
a mass limit may appear to be a long winded way of mak- 
ing this calculation, as it is straightforward to analytically 
calculate these numbers for a given ICMF. However ana- 
lytical calculations do not give an estimate of the error as- 
sociated with stochastic sampling of the ICMF, an effect 
that increases as you estimate the fraction of mass above 
an increasing mass limit. As estimating realistic errors for 
this and similar studies is one of the goals of this study, we 
choose to use the monte carlo technique in the subsequent 
analysis. 

With the ages and mass of clusters for NGC 3256 calcu- 
lated we found the total mass of clusters less than 10 million 
years old, and with a mass greater than 1O 4,7 M0. The lower 
mass limit of 1O 4 ' 7 M0 was chosen as a conservative esti- 
mate to ensure our sample was complete for this range of 
masses and ages. This gave a mass of 1.66 x 10 7 Mq con- 
tained within clusters. The fraction of mass expected above 
this mass limit was taken from our model cluster popula- 
tions to be 0.61 ± 0.09. We expect the total mass contained 
within all clusters in this age range tobe2.73±O.42xlO 7 M . 
This was divided by the age range of 7 million years to find 
the stellar mass forming in clusters per year, 3.9O±O.6OM0 
yr _1 . We used an age of 7 rather than 10 million years as 
extremely young clusters with ages less than 3 million years 
will still be embedded, an thus undetectable. 

We compare the mass of stars forming within clusters 
to the SFR which was calculated based on the IR luminos- 



10.0 r" 



ity from Sargent et al. 



46.17 Mq yr using the conversion of 



( 1989|) and con verted to a SFR o f 

(1998a I. 



The conversions given by Kennicutt ( 1998a I are based on a 



Kennicutt 



Salpeter IMF. In fitting SSP models to our photometric data 
we used SSP models which adopt a Kroupa IMF. This has 
the effect of underestimating the mass of clusters compared 
to a similar result using SSP models based on a Salpeter 
IMF. If we are to compare the mass of stars in clusters in 
clusters to the SFR both results must be made using the 
same IMF. We correct the masses of clusters for this by 
multiplying by 1.38. 

For NGC 3256, we find that the fraction of stars which 
remain within optically selected clusters after the embedded 
phase, F, is 0.12 ± 0.02. This is, however, a lower limit as 
various selection effects and fitting artefacts must be taken 
into account, which will be explored below. The implications 
of this result are discussed in § [5] 



3.2 Sources of Uncertainty 

In calculating V there are several sources of error we should 
consider. Here we focus on the errors associated with our 
study of NGC 3256 and also address issues which may arise 
with other data sets from different galaxies. 





Figure 6. Plot showing the effect of varying the mass function 
index on the fraction of mass above a mass limit (M > 10 4,7 ). 
The top panel shows the number of clusters found above the mass 
limit as a function of the ICMF index (—a), normalised to the 
value for a slope of a = 2.0. The middle panels displays R a , the 
ratio of mass found in clusters above the mass limit to the total 
mass of all clusters. The bottom panel shows R a normalised to 
i*2.0- 



3.2.1 Error Associated with the slope of the mass function 

As we previously mentioned, the slope of the cluster mass 
function was taken to be —2.0. However this may vary be- 
tween galaxies and has an effect on the fraction of mass 
we infer to be above the mass limit. With a steeper slope 
fewer high mass clusters are present so the fraction of mass 
above a mass limit is reduced. Conversely a shallower slope 
creates more high mass clusters and so the fraction of mass 
above a mass limit is increased. We have measured this effect 
by making sample cluster populations with a mass function 
slope between —1.8 and —2.2. The results of these simula- 
tions are shown in Figure [6] in which we assume a mass cut 
of M > 10 47 , the same used in our study of NGC 3256. 
The top panel in Figure [6] shows how the number of objects 
found above the mass cut varies with differing mass function 
slopes. 

The middle panel of Figure [6] shows the ratio of mass 
within clusters with masses greater than the mass limit com- 
pared to the total mass of all the clusters, R a . Over the range 
of 1.8 < a < 2.2 this ratio changes dramatically from 0.93 
to 0.25 respectively. Measured values of a vary but the best 
estimates are close to 2.0. The bottom panel in Figure [6] 
shows how the value of R a compares to the value for a mass 
function with an index of —2.0. If the cluster IMF index 
varies between 1.8 < a < 2.2 then we may be under or over 
estimating V by roughly ±50% by assuming a single value 
of a = 2.0. 
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Table 2. Parameters for the uncertainty of the magnitudes: 
Am„ = 10 d l+ d 2Xm,^ 



4.0 

log(M lim /M s ) 



Figure 7. The fraction of mass in clusters above a mass limit 
(Mjwn) f° r an ICMF based on a Schcchtcr function with /3 = 2.0 
and variable M* value. The black line shows the same fraction for 
a power law ICMF with identical slope, as well as three Schechter 
functions with differing M„ values (dashed line, 10 5,2 Mq; solid 



line. 



1O 6 O M ; dot-dashed line, 10 6 - 8 M Q ) 



3.2.2 Power Law or Schechter Function? 

The exact form of the initial cluster mass function (ICMF) 
is still under debate. In the case of NGC 3256 we have as- 
sumed a simple power law of the form shown in Equation [2] 
For other galaxies an index of a w 2 has been found (e.g. 
Zhang fc Fall||1999| |de Grijs et al.|[2003cl |Bik et al.|[2503j 
McCrady fc Graham|2007 1 and is expected from theoretical 
considerations (Elmcgreen & Efremov 19971. Recent stud- 
ies have, however, shown that there may be truncation in 
the mass function at the high mass end ( jGieles et alT20 06b 



[Bastian 2008 Larsen 2 009[). This can be represented by a 
Schechter function ( Schechter 1976 1 which behaves as a pure 



power law in the low mass regime and has an exponential 
fall off above a given value, M*. 

In Fig. [7] we show how the resulting fraction of mass in 
clusters varies for different mass cuts and different trunca- 
tion values, M*. Overall, the difference between the results 
with an underlying power-law or a Schechter function are 
small. In the specific case of NGC 3256, if a truncation ex- 
ists it is likely to be well above 10 6 Mq, hence any effect on 
the derived F values will be minimal. 



3.2.3 Error Associated with Fitting Cluster Parameters 
using SSP models 



Bastian et al. ( 2005a I and Gieles et al. ( 2005 \ compared the 



results of SSP fitting for clusters in M51 with the results 
for a model population, noting that the number of clusters 
at certain ages and masses can be enhanced or reduced by 
the fitting method. This can affect the calculated total mass 
of clusters when extrapolating from the mass found in high 
mass clusters. To investigate this we made a model popula- 
tion of clusters based on the solar metallicity GALEV SSP 
model tracks (Padova isochrones) . Ages of clusters were ran- 
domly distributed in the range < t < 10 10 yrs. The masses 
of clusters were assigned stochastically assuming a power law 
distribution of slope —2.0. Based on the age, masss and dis- 
tance to NGC 3256 we then assigned magnitudes from the 
SSP models. 

Photometric errors for each band were estimated by ex- 
amining the errors of the clusters measured in our sample 
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Figure 8. Density map showing the difference in cluster param- 
eters between our model cluster population and the results of the 
cluster fitting procedure. Positive values indicate values where we 
over estimate the number of clusters. The inset panel shows the 
mass of the input model population compared to the mass ob- 
tained by the cluster fitting procedure for various mass cuts for 
clusters with ages less than lOMyr. 



from NGC 3256. We were able to approximate the error with 
the relation Am„ = 10 dl+d2 Ym " . The values of di and d 2 are 
displayed in Table [2] A random correction to the magnitude 
in each waveband was added in the range -AM„ to +AM„. 
We produced a catalogue of 2607 clusters which fulfilled our 
magnitude limit. 

We ran the same fitting procedure on this artificial cat- 
alogue as we had done on the data for NGC 3256. Fitting 
cluster photometry to the solar metallicty Padova SSP mod- 
els. When examining the results of this fitting procedure we 
were as stringent with the artificial clusters as we were with 
the real data. We only accepted those results which met 
our criteria to be a good fit, as defined in § |2.3| We then 
compared the known mass, age and extinction of the model 
population to the results from the model fitting. We also 
excluded those model clusters with ages less than 3 million 
years, assuming that in reality these young clusters would 
be embedded and unobservable. 

We have compared the age-mass distribution of the 
model cluster population to the results of the SSP fitting 
routine in Figure [8] The colour scale represents the differ- 
ence in the number of clusters found with a given age and 
mass compared to the input model population. The white 



© 2009 RAS, MNRAS 000,[T]fT6] 



On the Fraction of Star Clusters Surviving the Embedded Phase 9 



slope across the Figure shows a large deficit of clusters, this 
is due to the stringent limits we apply to our clusters, accept- 
ing only clusters with low errors (< 0.2). Any faint clusters 
will have large errors and so we reject the fit. Figure [8] shows 
that young clusters are not fitted as we assume we will not 
detect young embedded clusters. The inset in Figure[8]shows 
a comparison of the mass in clusters above a mass limit, com- 
paring the results of cluster fitting to the known values of the 
input model cluster. For mass cuts below 10 5 Mq we tend 
to underestimate the mass above the mass cut. This implies 
that measuring the total mass of stars within clusters would 
be underestimated due to the effects of the fitting proce- 
dure. This effect is small however and for the mass cut used 
in our study of NGC 3256 results in an underestimation by 
5%. For the largest mass cut used we in fact over estimate 
this mass but only by 10%, though it should be noted with 
higher mass cuts we are more likely to encounter errors due 
to the stochastic sampling of massive young clusters. 

The test we have performed here to estimate the error 
associated with the SSP fitting does not take into account 
the intrinsic problems with either the Padova or Geneva SSP 
models. These synthetic stellar evolution tracks do not per- 
fectly represent real star clusters, especially at young ages. 
Even for models without stellar rotation or binary stars the 
uncertainties associated with massive stars can be large, and 
when these effects are included the errors increase. The true 
affect of fitting to inaccurate SSP models is ultimately hard 
to quantify and the analysis we have carried out can only 
be considered as a rough estimate of this effect. 
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Figure 9. The effect of under or over estimating the metallicity 
of clusters on the measurement of T. The red line represents the 
results of fitting a population of solar metallicity with twice solar 
metallicity SSP models, while the blue line shows the opposite 
(fitting a twice solar metallicity population with solar models). 
The green line is the same as the inset of Fig. [8] 



3.2.4 Error Due to an Unknown Metallicity 

In determining the cluster parameters of age, mass and ex- 
tinction we have fitted the photometric data to a SSP model 
of a certain metallicity. Using optical spectroscopy, T07a 
found that the metallicity of clusters and Hn regions in 
NGC 3256 was between 1 and 1.7 Z . The Padova SSP 
models used have only two metallicities in/near this range, 
namely solar and twice solar. Here we test the effect of our 
metallicity assumptions on our derived value of Y. 

Just as we did in § |3.2.3| we produced model cluster 
populations based on the solar metallicity Padova tracks, 
however we then tried to fit these data to the twice solar 
Padova SSP tracks. We also produced a model population 
based on the twice solar Padova tracks and fitted them to 
the solar Padova tracks in an attempt to examine the effects 
of over or under estimating the metallicity. These results are 
shown in Figure[9] the red line shows what happens when we 
overestimate the metallicity, using a twice solar metallicity 
track to fit clusters with a solar metallicity. In this case we 
overestimate the mass in cluster above a particular mass 
cut. This overestimation changes a little depending on the 
particular mass cut but is approximately a factor of two for 
the case of NGC 3256. 



3.2.5 Error Due to Selection Effects 

In measuring V we need to know how efficient we are at de- 
tecting massive, bright clusters, which we use to infer the 
total mass of stars in clusters. Measuring this efficiency is 
not trivial and a rigourous, quantitative discussion is not in- 
cluded here. Instead we examined our HST images of NGC 
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Figure 10. Colour composite image of a particularly crowded 
section of NGC 3256. Objects detected in the B band are circled 
in green. Objects in yellow boxes passed our criteria for a good fit 
to the SSP models. 



3256, overlaying the objects selected from the B band im- 
ages. Assuming our detection was perfect all the brightest 
objects should have been identified. In reality this is not pos- 
sible as dust will obscure some bright clusters and crowding 
makes identification of individual clusters difficult. 

We selected the SExtractor detection parameters to 
minimise the chance of two close clusters being identified 
as a single object. A typical field from a crowded region 
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Table 3. A list of the possible sources of uncertainty in measuring 
r for NGC 3625. The effect of each source is expressed as the 
factor applied to correct the measure value of T. 



Source 



Variation 



Effect 



CMF slope (a) 
Schechter Function 
SSP Fitting 
Metallicity 
Selection 



a ■ 



1.8 

: 10 5 



- 2.2 
- 7 Mr. 







O.5Z - 2Z C . 



0.66 - 2.38 
3.72 - 1.33 

1.04 
0.5 - 2 

1.25 



of NGC 3256 is shown in Figure |10| and we have overlaid 
the detected clusters with green circles, those clusters which 
passed our criteria for a good fit have been identified by yel- 
low boxes. The brightest of clusters are all selected and have 
good fits to SSP models. There are several clusters which we 
detect but do not have good fits, on examining these clusters 
we find that they have either one or two wavebands with er- 
rors greater than 0.2 mag and are thus rejected, and labelled 
as poor fits. 

We see from Figure [10] that the actual detection of clus- 
ters is not limiting our cluster sample, the ability to accu- 
rately fit SSP models is a larger effect. However in our syn- 
thetic cluster population models we apply the same criterion 
to determine good fits from poor fits. We are very severe in 
the limits we apply, disregarding any cluster which has an 
error greater than 0.2 mag in any of the 4 wave bands, ft is 
this criterion which limits our final sample, and is taken into 
account by modelling cluster populations. We can test this 
by looking at the fraction of bright, young, massive clusters 
we reject as poor fits in both the real sample and synthetic 
cluster populations. We examined the percentage of clusters 
with accepted good fits, for clusters younger than 10 Myr 
and with masses greater than 10 4 ' 7 Mq. In our sample from 
NGC 3256 this percentage was 47% whilst in our synthetic 
population the result was higher at 76%. 

The discrepancy in these results could have several 
causes, firstly our synthetic clusters have a more conser- 
vative estimate on the level of dust extinction than we 
might expect. With a greater level of dust extinction clusters 
become fainter, photometric errors become larger and the 
chance we reject the cluster fit increases. Secondly the syn- 
thetic cluster population is entirely theoretical and so poor 
photometry from overcrowding or from misaligned apertures 
is not taken into account. We can only estimate the percent- 
age of clusters we either miss altogether or discount due to 
a poor cluster fit. An estimate of 50% is the worst case sce- 
nario but some of this is accounted for by our model cluster 
populations to which we apply the same stringent criteria for 
good cluster fits. Taking this into account we estimate that 
roughly 20% of clusters are either not detected or rejected 
as poor cluster fits in addition to the fraction we account for 
based on the synthetic cluster population models. 

3.3 The Value of V in NGC 3256 

We are in a position to fully calculate the fraction of stars in 
clusters (r) for NGC 3256, including all the possible souces 
of error. Initial inspection of our data for NGC 3256 revealed 
a value of 0.117. We must account for the over or underes- 



timation introduced when we fit data to cluster models. In 
§ |3.2.3| we calculated that we would underestimate the mass 
above our mass limit by 5%, however to make this calcula- 
tion we made model cluster photometry based on SSP mod- 
els, and then compared this to the same models. In reality 
clusters are not perfectly modeled by any SSP and this un- 
derestimation is likely larger in reality. We have been rather 
sceptical in assuming that we only recover 80 ± 10% of the 
mass when we use this cluster fitting method. This increases 
the value of F from 0.12 ± 0.02 to 0.15 ± 0.03. 

If we assume all our clusters were solar metallicity the 
we expect to recover 96% of the mass using a mass cut of 
10 4 ' 7 Mq. A more likely scenario is that the cluster popula- 
tion has metallicities in the range 1—2 Zq, this will affect the 
total mass estimated to be in clusters compared to the ac- 
tual mass, as shown in Figure [9] As the metallicity of each 
individual cluster is unknown we cannot be certain if the 
mass in clusters calculated is an under or over estimation of 
the true mass. We expect clusters in NGC 3256 to be closer 
to a solar metallicity than twice solar, and so estimate that 
our cluster fitting method underestimates the total mass by 
a factor of 0.8±q;| taken from figure [9] for the mass cut used 
in our study. This further increases T from 0.15 ± 0.03 to 

o.i8±8:8l- 

We also need to consider the fraction of clusters that 
were not detected or thrown out as poor cluster fits. We 
estimated this fraction to be 20% in section [3.2.5| Combining 
these effects we arrive at a more assured value for T than the 
value the data initially suggests. The value used throughout 
the remainder of this study is 0.23±o;?q. 



4 RESULTS FROM OTHER DATA-SETS 

Fitting photometric observations to cluster evolutionary 
tracks is very time consuming. This study does not attempt 
to replicate the workings carried out on the clusters of NGC 
3256 on numerous other galaxies. Instead we have used pre- 
vious results from other studies, in an effort to calculate 
the ratio of stars forming within clusters. The various data 
sources are discussed below. 



4.1 NGC 1569 



We use the data set of Anders et al. (2004b I, which includes 



ages and masses for 161 clusters in the galaxy NGC 1569. 
Just as for NGC 3256 we consider clusters younger than 10 
Myrs. W ith the small distance to NGC 1569 1 Anders etUI, 
(2004b| was able to observe many low mass clusters, and we 
can use a low mass cut to estimate the total mass in clusters. 
In order to avoid measurements which are affected heavily 
by stochastic sampling of the ICMF we must use the lowest 
mass limit possible whilst being wary of the observational 
limits. We chose three lower mass limits, 10 30 M , 10 32 
Mq and 10 3 ' 4 Mq. In the same manner as we did for NGC 
3256 we extrapolated the mass found above this limit up to 
the total mass in clusters, this resulted in consistent answers 
from all three of these limits. Averageing these three results 
we calculated the total mass in clusters to be 3.52 ± 0.05 x 
10 5 Mq . Assuming we are unable to observe any embedded 
clusters younger than three million years the resulting CFR 
is 0.05 Mq yr _1 . The star formation rate is taken from an 
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Ha measurement from Moustakas & Kennicutt ( 2006 1 and 
the Ha SFR calibration in Kennicutt (1998a I, using the 



same adopted distance to NGC 1569 as Anders et al. ( 2004b ) 



of 2.2Mpc. Giving a SFR of 0.36±0.02 M yr~ 
the ratio of the CFR to the SFR, 13.9 ± 0.8% 



4.2 NGC 6946 



T is simply 



We used data from |Larsen| ( |2002[ ) who quoted ages and 
masses for 90 clusters in NGC 6946. Considering clusters 
younger than 10 Myrs and a lower mass limit of 10 3 ' 5 Mq 
we calculated the mass in clusters above this limit, giving 
1.284 x 10 5 Mq. The fraction of mass expected above this 
mass limit is calculated as in previous sections and results 
in a total inferred cluster mass of 1.95 ±0.23 x 10 5 M . We 
continue to assume that clusters younger than 3 Myrs are 
embedded an unobservable and so have an age range of 7 
Myrs, thus givi ng a CFR of .022 ± 0.003 M yr" 1 . 

The data of Larsen ( 2002 ) only comes from one pointing 



of the HST WFPC2 camera and does not include the whole 
galaxy. To calculate V we need the SFR across this region 
only. This was found using the area covered by the WFPC2 



and the SFR density taken from Larsen ( 2002 1, giving a SFR 
of 0.1725 M yr -1 . Consequently T is 12T5=h|;| %, although 
we note this HST pointing included the centre of NGC 6946 
and so the SFR might be higher that the global average we 
used. 



4.3 Small Magellanic Cloud 



We have used the data set generated by Hunter et al. ( 2003 1, 



a catalogue of ages and masses for 191 clusters in the SMC. 
This catalogue is thought to be incomplete for clusters 
younger than 10 Myrs, so we calculated the mass in clus- 
ters between 10 and 1 00 Myrs. With the S MC lying at a 
distance of 61 ± 3 kpc ( jHilditch et al.|[2005[ ) the data set is 
complete to very low cluster masses. 

Using lower mass limits of 10 26 , 10 2 8 , 10 30 & 10 32 M 
to calculate the mass in total mass in clusters gave consistent 
results of 1.59±0.03 x 10 5 M . We assumed a constant CFR 
and SFR over the age range 10-100 Myrs, giving a CFR of 
1.77 ± 0.03 x 10 -3 M yr" 1 . The SFR which we assume to 
be constant over the last 100 million years is taken from the 
extinct i on co rrected Ha luminosity value of |Kennicutt fc| 
Hodge| |1986[ ), converted to a SFR of 0.043 M yr -1 using 
the Equations given by Kennicutt ( |199 8a). The resulting 



value for F is 4.2 ±g;§ % 



Our value for F agrees extremely well with the value 



calculated by Gieles & Bastian (2008) of 3-5%, derived using 



size-of-sample effects on the same cluster sample. |Kruijsse"n| 
& Lamers ( 2008 1 looked at the ratio of clustered to field stars 



in the SMC and quote a minimum F of 0.5% and suggest a 
more reasonable value of 10% as their best result. 



4.4 Large Magellanic Cloud 

Ages and masses for 748 clusters in the LMC were taken 



from the data set of Hunter et al. ( 2003 1. Just as in the case 



of the SMC this catalogue is believed to be incomplete for 
clusters younger than 10 millions years and so we studied 
clusters in the age range 10 — 100 million years. 



With the LMC being a relatively nearby galaxy at a 
distance of only 58.5 kpc ( Macri et al. ||2006 1 the sample is 
complete to low mass and we calculated the mass of clusters 
above four mass limits 10 2 ' 6 , 10 28 , 10 3 & 10 3 ' 2 M . The 
fraction of mass expected above these limits was found as 
described in previous sections, resulting in a total inferred 
mass in clusters of 6.31 ± 0.20 x 10 5 M . We assume a con- 
stant SFR and CFR over the 90 million year period consid- 
ered, resulting in a calculted CFR of 7.01 ± 0.22 x 10 -3 M 

yr- 1 . 

The SFR (0.12 M yr" 1 ) and the galaxy area (79 kpc 2 ) 
are both taken fromlLarsen & Richtler|(|2000|). F for the LMC 



was thus found to be 5.8 ± 0.5%. 



4.5 The Milky Way (The Solar Neighbourhood) 



Lada & Lada (2003) estimate the rate of embedded clus- 



ter formation within 2.0kpc of the sun to be between 2 — 4 
Myrs -1 kpc -2 , and claim the average mass of an embedded 
cluster to be 5OOM . Lada & Lada (20031 go on to esti- 
mate the number of embedded clusters that survive, with 
approximately 7% surviving to the age of the Pleiades. 

We convert the rate of embedded cluster formation to 
a mass formation rate assuming the average cluster mass of 
5OOM , giving O.15M yr -1 in the solar neighbourhood. We 
calculate a value for F by assuming the SFR is equal to the 
rate of embedded cluster formation. Comparing this to the 
mass found in older clusters, 7% though estimates of this 
percentage varies between 4% and 14% depending on the 



numbers used. Roberts (19571, Miller & Scalo (1978), and 



Adams & Myers (2001 \ also derived a similar percentage of 



star-formation in clusters for the solar neighbourhood. 



4.6 M83 



|Harris et al.| ( |2001[ ) presented a catalogue of 45 massive star 
clusters in the centre of M83, with reliable ages, masses and 
extinction values. Using this catalogue and a lower mass 
limit of 10 2 - 8 M we calculated the total mass in clusters 
younger than 10 million years to be 7.24 ± 0.23 x 10 5 Mq. 
To calculate the CFR we assume that clusters younger than 
3 million years old are embedded and so are unobservable. 
The resulting CF R is 0.10 ±0.01 M yr -1 . 

The study of |Harris et~aTl p00T| | only covers the cen- 
tral region of M83 and so we used an Ha image of M83 
to measure the Ha from the same region studied by |Har-| 
|ris et aL| ( |2001[ ) . We used archival narrowband and R-band 



images of M83 taken as part of the SINGS survey (Meurer 



et al. 20061. We measured the background subtracted flux 
for an identical area as that studied by |Harris et al.| ( |2001[ ) 
and converted this flux to a SFR using the calibrations of 



Kennicutt (1998a I, resulting in a SFR of 0.23 ± 0.03 M 



yr . We correct this value for extinction using the redden- 



ing curves of Calzetti ( 2001 1. The extinction value was taken 



as the averaged extinction of all clusters in the |Harris et al.| 
( |2001[ ) sample, giving a correction of 1.71, the resulting dust 
corrected SFR is 0.39 ± 0.06 M yr -1 . 

With the CFR and SFR known we calculate the value 
of T to be 26.7 ±|;g %. 
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4.7 The Antennae Galaxies 

Ages and masses have been published for 752 clusters in the 



Table 4. Summary of results included in this paper. 



Antennae Galaxies by Anders et al. (20071. We attempted to 



apply the same calculations as we did for other cluster pop- 
ulations, as described above. However using different lower 
mass limits we obtained a range of values for F, from 60% 
to > 100%; assuming a star formation rate of 20 M yr 
( Zhang et al.|206T l. Obviously values in-excess of 100% are 
unphysical and this problem is exasperated by the range of 
quoted SFR's for the Antennae Galaxies which can range 



from 5-20 M yr _i (Zhang et al 



2003) 



2001 



Knierman et al. 



We investigated the slope of the mass function for young 
clusters and found that this data set presents a very shallow 
slope of a ~ 1.6. Given that mass function slopes are usually 
in the range 1.8 < a < 2.2 (de Grijs et al. 2003a I a slope 



this shallow represents a very unusual population, unfortu- 



nately Anders et al. (20071 do not calculate a similar mass 
function for these clusters. Correcting for this shallow mass 
function gave results for T that were still above 50% and 
did vary significantly with the assumed lower mass limit. In 



comparison we noted that Fall et al. (2005) conclude that 
at least 20% but possible all clusters form in clusters; the 
value for F may well be high for the Antennae. 

Given the range of values we can calculate for F based 



on the Anders et al. ( 2007 I data set and the extremely shal- 
low mass function of young clusters in this sample combined 
with the range in quoted SFR's for the Antennae we do not 
publish a value for this system. It serves as a reminder of 
how difficult these calculations can be and how dependent 
they are on accurate knowledge of the SFR. 



5 THE VARIATION OF V WITH SFR 

In Table[4]we present all the information used hereafter, in- 
cluding the values of F, star formation rate (SFR) and area 
over which these quantities were measured. The superscript 
p after the galaxy names indicates that only a partial area 
of the galaxy was used to derive these results, as described 
in the previous sections. At this point it is worth noting that 
because results have been obtained from various data sets 
the methodology used differs between results. Overall the 
consistency of results is thought to be robust as measure- 
ments of r are based on the brightest and easiest clusters to 
detect which should give consistent results. It is harder to 
estimate the errors associated with F based on other data 
sets. Results shown in Figure[TT]for galaxies other than NGC 
3256 and the Milky Way only include errors associated with 
the uncertainty in the fraction of mass expected above what- 
ever lower mass limit was used. Errors due to uncertainties 
in the metallicity and SSP fitting have not been estimated 
as this is difficult without reprocessing the entire data set. 

Figure |TT| shows the relation between the SFR density 
(Esfr) and F for all galaxies discussed in this paper. In ad- 



dition we plot results for 3 galaxies from Gieles (2009a I in 
which CFRs (and subsequently V) are calculated via com- 
parison with empirical luminosity functions. It is immedi- 
ately obvious that V increases with the SFR density. Over 
three orders of magnitude a power-law relationship holds 
and the dot-dashed line shows a least-squares fit to the data 
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Figure 11. Graph showing V against star formation rate density 
for all galaxies examined in our study. The dot-dashed line shows 
a least-squares best fit power law to the distribution. Diamonds 
show data points taken from |Gie lcs (2009a) and are not included 
in calculating the best fit to the data. 



of the form, Y oc Y>g FR . The numerical version of the derived 
relationship is displayed in Equation [3] 

T(%) = (29.0 ± 6.O)E£ 24 ± o ' o4 (M yr- 1 kpc~ 2 ) (3) 

From this Equation we can predict at what density we 
might expect to see all stars in clusters, effectively a T of 
100%, this would occur at a SFR density of 7.5 x 10 3 M 
yr -1 kpc -2 , a rather high value. Interestingly we can in- 
tegrate this value over the time it takes a cluster to form, 
assuming this takes roughly three million years we expect 
2.25 x 10 4 M pc -2 , which is approximately the density of 
a typical cluster. 



Although the results taken from Gieles ( 2009a I are not 
calculated in the same manner as we have for other data sets 
these points fit the trend quite well, particularly at low SFR 
densities. It should be noted that our calculation of the best 
fit power law does not include these additional data points, 
given the different method used to determine these results, 
they are shown here to demonstrate the consistency of our 
results. 

We have investigated whether the observed correlation 
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(Eqn. [3| could be due to a systematic effect, namely if the 
mass function index varied with environment. Assuming the 
LMC cluster distribution is well approximated with a — 2, 
an index of a ~ 1.2 would be required in order to bring the 
results of NGC 3256 into agreement with the LMC. This 
kind of gross deviation is ruled out by our direct observations 
of the cluster population (e.g. Fig. [4J. 

A similar conclusion (that more stars are formed in 
bound clusters with increasing SFR density) has been 
reached based on the fraction of light observed in clusters 
relative to the host galaxy. Specifically, Meurer et al. (1995) 
and Zepf et al. (1999) found this for merging/starburst 
galaxies and Larsen & Richtler ( 2000 1 and Larsen ( 2004 1 



found this for a larger sample that includes starburst and 
quiescent galaxies. While measuring the fraction of light in 
clusters is significantly easier than deriving V, hence suit- 
able for construction of large surveys, it is affected by the 
star formation history of the galaxy and possibly by differ- 
ential extinction effects. The exact relation between F and 
the fraction of light observed in clusters will be modelled in 
a future paper. 

Kennicutt ( 1998b I investigated a global Schmidt- 



Kennicutt law across galaxies, correlating the SFR density 
to the surface gas density. As V appears to follow a power- 
law relationship with the SFR density using the Schmidt- 



Kennicutt law defined in Kennicutt ( 1998b I we can con- 



struct a relation between F and the surface gas density, 
shown in Equation [4] 



r(%) = (4.1±1.9)E° a 3 s ' 



7 {M Q pc 



(4) 



We show the derived surface gas density corresponding 
to a given SFR density on the top axis of Figure |11| Al- 
though our data only charts the value of F for gas densities 
in the range 1 — 300 Mq pc -2 it is insightfull to interpret 
this relation to the more extreme star forming environments. 
Arp220 represents one of the most active star forming galax- 
ies found \ Wilson et al.|2006[ ), with a global SFR of 240 My 
yr _1 based on the FIR luminosity of panders et al. (20031 
and the calibration of |Kennicutt| ( |1998a[ ). Assuming that 



Arp220 occupie s an area of roughly 1 kpc ( Anantharama- 



iah et al. 



2000), and thus has a surface area of 1 kpc on 



the sky we can derive the SFR density is approximately 240 
Mq yr _1 kpc -2 . Such a high SFR density translates to a 
value of r of roughly 85% using Equation [3] indicating that 
almost all the star clusters in Arp220 would be expected to 
survive the embedded phase, and the CMF would be almost 
identical to the true underlying initial CMF. Hence, there 
are regions where we would expect most/all young stars to 
be found in clusters. 



how r varies with the star formation rate of the host galaxy 
and we summarise our conclusions below. 

(i) The cluster formation history of NGC 3256 shows an 
increase in the star formation rate over the last 100 million 
years, most likely caused by the ongoing merger of the two 
progenitor spiral galaxies. The CMF for NGC 3256 is best 
described by a power law with slope a = 1.85 ± 0.12. 

(ii) r may be calculated directly if an accurate cluster 
population is known, however there are several possible 
sources of error in making these calculations. The effect of 
SSP fitting must be taken into account and an unknown 
metallicity may produce additional uncertainties. Selection 
effects may also alter the value of F but this can be esti- 
mated. The calculated F does depend on the form of the 
cluster mass function used, we assume a simple power law 
function with slope a — 2.0. We have shown that a Schechter 
function produces similar but higher values, however this ef- 
fect is small if the M*is high (M*> 1O 6 M ). 

(iii) We found a weak positive correlation of F with the 
total star formation rate. We find a strong relation between 
F and the star formation rate density, which we have written 
as a power law type relation, F(%) = (25.5 ± 6.0)S° 4 2 R ±0 05 



Larsen 



(Mq yr -1 kpc --3 ). This is similar to that found by 
|fc Ric htler ( 2000| for the fraction of U-band light in clusters 
relative to the total galaxy. This result can also be inter- 
preted as a correlation with the surface gas density through 
the Schmidt-Kennicutt Law ( |Kennicutt| 1998b| . This implies 
that either clusters born in high density environments are 
more resistant to disruption in the embedded phase, or the 
environment changes the fraction of stars born in clusters. 

Measuring F is not trivial and it does have several 
sources of error which must be taken into account. How- 
ever it is possible to make these calculations which may be 
further refined with consistent measurements using the same 
instruments, detection methods and cluster fitting models. 
This study is intended to show that if all sources of uncer- 
tainty are taken into account and measured then calculating 
r is possible. Future work, including larger homogeneous 
samples and deeper observations to constrain the form of 
the cluster mass function to lower masses, will be needed 
to conclusively investigate the intriguing relation between F 
and the SFR density of the host galaxy. 
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6 CONCLUSIONS 

Over the course of this paper we have shown that it is 
possible to accurately measure the fraction of stars found 
within young clusters, a parameter we have termed F. This 
is achieved by obtaining ages and masses for the cluster pop- 
ulation of NGC 3256 through multiple waveband photome- 
try and comparison with synthetic stellar population mod- 
els. We have also been able to calculate F for several other 
galaxies using published cluster populations. We examined 
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